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Abstract 

A frequently desirable characteristic of chemical kinetics systems 
is that of persistence, the property that if all the species are initially 
present then none of them may tend toward extinction. It is known 
that solutions of deterministically modelled mass-action systems may 
only approach portions of the boundary of the positive orthant which 
correspond to semi- locking sets (alternatively called siphons). Conse- 
quently, most recent work on persistence of these systems has been 
focused on these sets. 

In this paper, we focus on a result which states that, for a conserva- 
tive mass-action system, persistence holds if every critical semi-locking 
set is dynamically non-emptiablc and the system contains no nested 
locking sets. We will generalize this result by introducing the notion 
of a weakly dynamically non-cmptiable semi-locking set and making 
novel use of the well-known Farkas' Lemma. We will also connect 
this result to known results regarding complex balanced systems and 
systems with facets. 
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1 Introduction 



An elementary chemical reaction is given by a set of reactant species react- 
ing at a prescribed rate to form a set of product species. Under appropriate 
assumptions (well-mixing, constant external conditions, large number of re- 
acting molecules, etc.) such systems can be modelled deterministically by 
an autonomous set of ordinary differential equations over continuous vari- 
ables representing the concentrations of the reactant species. The resulting 
mathematical models have a long history and enjoy applications in fields 
such as systems biology, industrial chemistry, atmospherics, etc. [UEHSIIII] 

One topic which has gained significant attention recently has been that 
of persistence [TH3ll5l [T^l[T6] . A chemical kinetics system is said to be persis- 
tent if no initially present species may tend toward extinction. In general, 
determining whether a system is persistent can be difficult, but it has been 
significantly simplified by recent work. In particular, in [3] the authors 
show that the boundary of the positive orthant can be divided into subsets 
Li (roughly faces of M>q) suc h that trajectories may only approach Lj if 
I C. {1, ... ,m} is a siphon (called semi-locking sets in this paper). Conse- 
quently, in order to determine persistence it is sufficient to consider only the 
behaviour near sets Lj corresponding to sets / which are semi-locking sets. 

In the same paper, the authors introduce a condition on the sets / called 
dynamical non-emptiability. This condition corresponds to a dominance 
ordering of the reactant monomials which compose the differential equations. 
Together with some technical assumptions, they prove that the dynamical 
non-emptiability of every critical semi-locking set is sufficient to guarantee 
persistence. While certainly powerful, the result is limited by its restrictive 
assumptions, which include that the system be conservative and have no 
nested critical deadlocks (Theorem 4, [3]). In this paper, we extend this 
result by introducing the notion of weak dynamical non-emptiability. Our 
methodology differs significantly from that used in [3] and will allow us to 
relax several assumptions while maintaining persistence. Our methodology 
depends crucially on a novel use of the famous Farkas' Lemma [6]. 

For many systems, proving persistence is tantamount to proving the 
global asymptotic stability of some positive equilibrium state, a character- 
istic which is typically highly desirable. One such class of systems are the 
complex balanced systems first considered in 7. !). 1 1 In [11] . the authors 
proved that relative to each stoichiometric compatibility class, the invari- 
ant spaces in which solutions are restricted, there is precisely one positive 
equilibrium state and that this state is locally asymptotically stable. 

Despite significant empirical support, however, the hypothesis that this 
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equilibrium state is in fact globally asymptotically stable relative to its com- 
patibility class remains unproven. In the literature, this hypothesis has been 
termed the Global Attractor Conjecture and a significant amount of work has 
been done attempting to confirm it [UE1EJCE31CL6] • Of particular interest to 
us is the work of [TB], where the authors prove that the w-limit set consists 
of either this unique positive equilibrium state or of some complex balanced 
equilibria on the boundary of the positive orthant. Consequently, for com- 
plex balanced systems, persistence is a sufficient condition for the global 
asymptotic stability of this equilibrium state relative to its stoichiometric 
compatibility class. It follows that the work of this paper affirms the Global 
Attractor Conjecture for chemical kinetics systems satisfying the relevant 
assumptions which are also complex balanced. 

Persistence is also considered in [2] where the authors prove that mass- 
action systems are persistent if every semi-locking set I corresponds to a 
facet of the compatibility class (i.e. a facet is a codimension one face of the 
compatibility class). The authors connect their work with [3] by showing 
that sets I corresponding to facets are dynamically non-emptiable under 
some assumptions of the system. We will prove that there is no restriction 
for weak dynamical non-emptiability; that, in fact, every set / corresponding 
to a facet is weakly dynamically non-emptiable. 

It should be noted that our notation and terminology will often dif- 
fer from that of [3]. We have attempted, wherever possible, to adopt ter- 
minology specific to chemical kinetics systems (for example, we follow pQ 
in naming siphons according to their chemical kinetics equivalent of semi- 
locking sets). Throughout this paper, we will let M> and M> denote the 
m-dimensional spaces with all coordinates strictly positive and non-negative, 
respectively. 

2 Background 

Within the mathematical literature, several distinct ways to represent chem- 
ical reactions have been proposed. In this paper, we will follow closely the 
reaction-oriented formulation of [3]. (For examples of species- and complex- 
oriented formulations, see [T7] and [11], respectively.) 

We will let Aj denote the species or reactants of the system and let S 
denote the set of distinct species of the system. We define |<S| = m. The 
reactions of the system will be represented as 

m m 

K t : ^ (1) 

i=i j=i 
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where a^, fyj G Z>o are called the stoichiometric coefficients and ki > is 
the rate constant associated with the i th reaction. We will let 1Z denote the 
set of reactions of the system and define \7Z\ = r. 

2.1 Mass- Action Kinetics 

It is convenient to collect the net gain or loss of each species as a result of 
a specified reaction in a matrix T G Z mxr where the entries of T are defined 
according to 

[T]ji = f3ij - aij (2) 

for i = 1, . . . , r and j = 1, . . . , m (notice the reversal of indices). The entries 
[r]jj can be intuitively interpreted as the change in the j th species as a result 
of each instance of the i th reaction. 

We will let Xj = [Aj] denote the concentration of the j th species and 
denote by x = [x\ X2 ■ ■ ■ x m ] T the concentration vector. We define the 
reaction vector -R(x) G M> according to 

m 
3=1 

for i = 1, . . . ,r. The terms i?i(x) correspond to the reaction rates of chem- 
ical reactions under the assumption of mass-action kinetics. (For further 
discussion on alternatives to mass-action kinetics, see [3].) 

Since represents the net stoichiometric change in the j th species as 
a result of the i th reaction and Ri (x) represents the rate of occurrence of the 
jth reac ti OI1) ^ follows that [r]jji?j(x) is the rate of change of the j th species 
as a result of the i th reaction. It follows that (prj is governed by the system 
of differential equations 

^ = f (x) = rfl(x). (4) 

It follows from (j3]) that solutions are not able to wander around freely in 
M m since f(x) G range(r) for all x G K> - Trajectories are in fact restricted 
to stoichiometric compatibility classes (see [llj). 

Definition 2.1. The stoichiometric subspace for a chemical reaction 
mechanism (CP is the linear subspace S C R m given by 

S = range(r). 

The dimension of the stoichiometric subspace will be denoted by s = rank(r) . 



4 



Definition 2.2. The positive stoichiometric compatibility class con- 
taining the initial concentration xq € M™q is the set C Xo = (xo + S) n M> . 

Proposition 2.1 ( pTT]fTTj ). Let x(i) 6e i/ie solution to ^ with x(0) = xq € 



So- 



T/ien x(t) € C Xo for t > 0. 



A significant amount of literature exists analysing and restricting the 

behaviour of solutions of @ ^000111111131121. We wil1 omit most of 
this discussion, instead focusing on the analysis conducted in [3] and the 
related paper pp. 

2.2 Conservative Systems and Semi-Locking Sets 

In this section, we introduce several background concepts from [3] which 
will be important to us. We have adapted some of their terminology to the 
more specific framework of chemical kinetics. 
We start by defining a conservative system. 

Definition 2.3. A chemical kinetics system is said to be conservative if 

there exists a c £ R> such that 

c T r = o T . (5) 

We will call any vector c £ R™ satisfying a conservation vector. 

An important property of conservative systems is that solutions of (jU) 
with xo £ M> are bounded and do not approach the origin. This follows 
from (j4j) and (|5j) since these collectively imply i [c r x(t)] = from which 
it follows that c T x(T) = c T x(0) > for all T > 0. Divergence of any 
species to infinity, however, implies the existence of a sequence {t n } such 
that c T x(i n ) — y oo as n — y oo, while convergence to the origin implies a 
sequence such that c T x(t n ) — >• as n — > oo. 

Related to the concept of a conservation vector is the concept of a semi- 
conservation vector. (These are called P-semiflows in [3] . We introduce the 
term semi-conservation to emphasize the connection with chemical kinetics.) 

Definition 2.4. Consider a nonempty set I C {1,2, ... ,m}. A chemical 
kinetics system is said to be conservative with respect to I if there exists 
a c € M>q satisfying 

Ci > 0, i € I , . 

Ci = 0, i#I {0> 
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so that 

c T r = o T . (7) 

We will call any vector c G M> satisfying |]5]) and ([Tj) a semi-conservation 
vector. 

Many of the results of this paper will depend on restricting our attention 
to only those components X\ of x G R m such that i E I for some set I C 
{1, . . . , m}. When / corresponds to only those indices for which ij / we 
will call I the support of x. 

The following concepts are the central focus of the papers pQ and [3]. 
They will be key to the analysis conducted in the rest of this paper. 

Definition 2.5. The nonempty index set I C {1, 2, . . . , m} is called a semi- 
locking set (alternatively, a siphon) if, for every reaction where an element 
Ai, i € /, appears in the product complex, an element from Aj, j € / 
appears in the reactant complex. A semi-locking set I is called a locking 
set (alternatively, a deadlock) if every reaction has an element Ai, i G I , 
which appears in the reactant complex. 

It is worth noting that the element Ai appearing in the product may 
be the same as the element Aj appearing in the reactant. The concepts 
of siphons and deadlocks are based on the theory of Petri Nets introduced 
in [12] and used extensively in [3]. We will, however, use the chemical 
reaction specific terminology of semi- locking and locking sets used in [1]. 

There are several intuitive consequences of semi-locking and locking sets. 
Since any species on the reactant side of a reaction being absent (i.e. Xi = 
for this species) implies by ([3]) that the rate of the reaction is zero, if all the 
species in a semi-locking set are absent then no species in the set will ever 
be produced (i.e. Xi(t) = for all t > for all species in the set). In this 
sense, the set of species is locked in place. In an analogous way, for a locking 
set all of the species are locked in place. 

The authors of [3] further classify semi-locking sets according to the 
following definition. 

Definition 2.6. A semi-locking set I C {1, 2, . . . , m} will be called critical 
if it does not contain the support of a semi-conservation vector. 

2.3 Persistence and Faces 

The primary goal of [3] is guaranteeing that solutions of do not approach 
cMIJ'q for any xq G M>q- I n other words, the authors sought to guarantee 



that if every species is initially present then none of them could not tend 
toward extinction as a result of the reaction system. 
The concept is succinctly stated as follows. 

Definition 2.7. A chemical reaction network with bounded trajectories is 
said to be persistent if, for any xo E R>o> it follows that w(xo) ncM> = 0. 

Persistence has historical roots in predator-prey modeling in biology but 
has been frequently applied to chemical reaction mechanisms in recent years. 
We require trajectories be bounded to avoid ambiguous limiting behaviour. 

In order to make the concept of persistence more manageable, we divide 
the boundaries of K> and C xo into faces. For technical reasons we will only 
be interested in the relative interior of these faces, which we will define as 
follows. (These sets are defined similarly in p], [3], and [5]. In [2], Lj is 
denoted Z[.) 

Definition 2.8. Given a nonempty index set I C {1,2, . . . ,m}, we define 
the sets Li and Fj to be 

Lj = {x E M> | Xi = for i E / and Xj > for i g" /} 
Fj = jx E C xo | Xi = for i £ I and X{ > for i g" /} 

The set Fj can also be given as Fi = C xo n Lj or Fi = (xo + S) D Lj. 

We notice that each x E cM^ can be placed into exactly one set Lj 
and each x E <9C X0 can be placed into exactly one set Fj so that these 
sets uniquely and completely decompose cM™ and <9C X0 , respectively. Not 
every /, however, necessarily corresponds to a Fj which contains points, i.e. 
it is possible that Fj = 0. Such sets will be said to be stoichiometrically 
unattainable. 

The following result corresponds to Proposition 1 of [3] and Theorem 2.5 
of pQ. This result places strong limitations on where w-limit points may lie 
on dM™ Q . 

Lemma 2.1. Consider a mass-action system and a nonempty index set 
I C {1, 2, ... , m}. If there exists a xo E M™ such that w(xo) n Lj ^ 0, then 
I is a semi-locking set. 

An important consequence of Lemma 12.11 is that, in order to prove the 
persistence of a mass-action system, it is now sufficient to prove only that 
uj(xq) n Lj = for all sets Lj corresponding to semi- locking sets /. In 
practice, this is a significant simplification. 
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2.4 Dynamical Non-Emptiability 

In Section 12.31 we saw that, in order to determine persistence of chemical 
kinetics systems, it is sufficient to consider behaviour near the sets Lj corre- 
sponding to semi-locking sets. To this end, Angeli et al. introduced several 
new concepts, including a notion of partial ordering on reaction rates and 
two cones: the feasibility cone and the criticality cone. These ideas culmi- 
nate in the concept of a semi- locking set being dynamically non-emptiable. 
Everything in this section can be found in [3]. 

We start by defining a partial ordering condition on the reactions of a 
system. 

Definition 2.9. Consider the nonempty index set I C {1,2, .. . ,m}. For 
lZi,lZj S 7Z, we will say that TZi =3/ TZj if > ajk for all k € I and the 
inequality is strict for at least one k € I. 

Intuitively, the partial ordering condition given in Definition I2.9I gives us 
an estimate on the magnitudes of the reaction terms Ri(x) near a set Lj. 
This is made explicit by the following result. 

Lemma 2.2 (Lemma 4, [3]). Consider a mass-action system and let I C 
{1, 2, . . . , in} be a semi-locking set. Suppose that TZi TZj- Then, for every 
e > 0, and each compact subset K of Lj, there exists a neighbourhood U of 
K in IR™ such that i?j(x) < ei?j(x) for all x G U . 

The concept of dynamical non-emptiability depends on two cones, the 
feasibility cone and criticality cone, which are defined as follows. 

Definition 2.10. The feasibility cone is defined to be 

T e (I) = {v G M> I Vi < evj, V Ki,Tlj G 1Z such that TZ { ^ TZj} 
where e > 0. 

Definition 2.11. The criticality cone is defined to be 

C(/) = {vei r > | [rv] fe <o, Vfce/}. 

We can now define dynamical non-emptiability, which is one of the major 
concepts introduced in [3]. 

Definition 2.12. A critical semi-locking set is said to be dynamically 
non-emptiable if there exists an e > such that 

T e (i)nc(i) = {o}. 
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2.5 Results of Angeli, De Leenheer, and Sontag 

The following two persistence results are proved in [3]. They are the basis 
of the original work contained in Section O (Two sets I\ and I2 are nested 
if h C I 2 or I 2 C h.) 

Theorem 2.1 (Theorem 2, [3]). Consider a chemical reaction network sat- 
isfying the following assumptions: 

1. the system is conservative; 

2. every semi-locking set contains the support of a semi-conservation vec- 
tor. 

Then the system is persistent. 

Theorem 2.2 (Theorem 4, [3]). Consider a conservative mass-action sys- 
tem satisfying the following assumptions: 

1. all of its critical semi-locking sets are dynamically non-emptiable; 

2. there are no nested distinct critical locking sets. 
Then the system is persistent. 

3 Original Results 

In this section, we generalize the results contained in Section 12.51 We start 
by presenting some necessary background material not presented in [3]. 

The following result is a version of the well-known Farkas' Lemma and 
should be contrasted with Lemma 5 of [3] as it will be used in similar fashion. 
(This formulation of the result follows from the statement of Farkas' Lemma 
given in [13] , taking a € W£ \ {0}.) 

Lemma 3.1 (Farkas' Lemma, [6]). Consider A € M mxn . Then exactly one 
of the following two conditions is true: 

1. There exists x S ]R> > x kei(A), such that Ax < 0. 

2. There exists y G M™ such that A T y > 0. 

The following persistence result can be found in [15]. We do not prove 
it here. 
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Theorem 3.1 (Theorem 3.13, |15j). Consider a general mass-action system 
with bounded solutions. Suppose that for every Lj corresponding to a semi- 
locking set I there exists an a satisfying 



and the following property: for every compact subset K of Lj, there exists a 
neighbourhood U of K in ]R> such that 



Then the system is persistent. 

3.1 Weak Dynamical Non-Emptiability 

In this section, we extend the notion of dynamical non-emptiability by mod- 
ifying the feasibility cone and introducing a kernel condition. We also refor- 
mulate the conditions required for inclusion in the feasibility and criticality 
cones as matrix conditions which will allow us to prove the main result of 
the paper (Theorem 13.51) . 

Our notion of dynamical non-emptiability depends on the selection of a 
set J C TZj where 



The key modification here is that we do not necessarily need to consider all 
pairs (i,j) satisfying TZi Ai TZj] often it is sufficient to consider a strict subset 
of these pairs of reactions. In Section [U we will see that this modification 
allow us to encompass more chemical kinetics systems than we would be 
able to otherwise. 

We will need the following two concepts. 

Definition 3.1. We define the feasibility cone relative to J to be 




(8) 



(a, x) < for all x € U. 



0) 



Tlj = {{i,j) € {l,...,r} x {l,...,r} | TZ t Aj Kj} . 



(10) 



T e (J) = {v G | Vi < evj, for all (i,j) € J} 



where e > 0. 



Definition 3.2. We define the kernel of I and J to be 



ker(J, J, e) = {v € R r > | [Tv] k = 0, for all k e I 

and Vi = evj, for all (i,j) £ J} 



where e > 0. 
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The following notion of dynamical non-emptiability is our own. It is 
more general than that contained in [3] in that it makes use of the freedom 
to select an appropriate J C IZi and broadens the inclusion principle to a 
kernel condition. (It is clear that the standard notion of dynamical non- 
emptiability is included as a special case of the following by taking J = TZj 
and recognizing that {0} C ker(/, J, e).) 

Definition 3.3. A critical semi-locking set I is said to be weakly dynam- 
ically non-emptiable if there exists an e > and a J satisfying &1(J\) such 
that 



In order to relate the above conditions to Farkas' Lemma (Theorem 13. 1|) 
we restate them as matrix conditions. We let nj = \I\ and nj = \J\. We 
define Ti £ I^ ?1 ^ XT ' to be the matrix T with the rows T^., k £ I, removed. 
We define Tj £ M. njXr to be the matrix where each row corresponds to a 
specific condition IZi =<!/ IZj, (i,j) £ J, so that in that row there is a one 
in the i th column, a — e in the j th column, and zeroes elsewhere. Lastly, we 
define f £ R( n i+ n j) xr to be 




The following result can be trivially seen. 

Lemma 3.2. The condition C(I) n J e (^) C ker(I,J,e) is satisfied if and 
only if 



3.2 Facets and Non-Critical Semi-Locking Sets 

In general, determining whether a semi-locking set / is weakly dynamically 
non-emptiable can be tedious. To this end, in this section we show that there 
are classes of semi-locking sets which are necessarily weakly dynamically 
nonemptiable: namely, semi-locking sets corresponding to facets (i.e. sets 
Fj of dimension s — 1) of a weakly reversible mechanism, and semi-locking 
sets which are non-critical. 

Facets are the central topic of consideration in [2], where the authors 
prove the following result. 

Theorem 3.2 (Theorem 3.4, [2]). Consider a weakly reversible mass-action 
system with bounded trajectories. Suppose that every semi-locking set I is 
such that F[ is a facet or empty. Then the system is persistent. 



C(/)nJ £ (J) C ker(/, J,e). 




>o 
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The authors also connect the notion of a facet with the traditional notion 
of dynamical non-emptiability (Corollary 3.5, [2]). Their result, however, 
overstates the implications of Fj being a facet. It can be shown that semi- 
locking sets I corresponding to facets Fj may fail to be dynamically non- 
emptiable if there is a reaction TZi such that: (1) IZi ^.j TZj for all IZj in the 
same linkage class (a connected portion of the reaction graph); and (2) IZi 
produces no stoichiometric change in the species in /. (See Example 1 of 
Section H) 

We now generalize this result by showing that there is no such exemption 
for weak dynamical non-emptiability. 

Theorem 3.3. Consider a weakly reversible mass-action system with a 
semi-locking set I C {1, . . . , m}. If Fj is a facet then I is weakly dynamically 
non-emptiable. 

Proof. We will follow closely the proofs of Theorem 3.2 and Corollary 3.5 
contained in [2]. 

In their proof for Theorem 3.2, the authors show that there exist Zj > 0, 
j G /, and ji G R, i G {1, . . . , r}, such that 



for all j G / and i G {1, . . . , r}. In other words, relative to the support of 
/, every reaction vector (the columns of T) lies within the span of a single 
vector which is strictly positive on the support of /. We will let z G 
denote the vector with the elements zj, j G /, indexed in order. 

It follows immediately from (jlip that every reaction in the system con- 
tributes either: (1) a net gain to all species in /, (2) a net loss to all species 
in /, or (3) no stoichiometric change to species in /. We can also divide 
the reactions according to the linkage classes k = 1, ...,£. By weak 
reversibility each linkage class is strongly connected. We will let de- 
note the reactions in the k th linkage class, and R+\ R_ , and Rq denote 
respectively the reactions in the k th linkage class which contribute a net 
gain, a net loss, or no change to all species in /. We notice that ji > for 
i G 7i < for i G R^ , and ji = for i G R { Q k) . Combined with (fTT]l . 
this division gives 



liZj — Pij 



(11) 




(12) 
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We now proceed to construct our set J C TZj according to (fTU|) . Since 
the system is weakly reversible, it follows that the product complex of every 
reaction TZi is itself a reactant complex for some other reaction (which we 
will denote TZi') in the k th linkage class. Since reactions may only produce 
simultaneous gain or loss to the species in /, it follows that: (1) TZ^ TZi 
if i G (2) TZi fa IZv if i G R^\ and (3) TZ V fa TZi and TZi fa TZ V if 

i G i?^- (Notice here that TZ V fa TZi and TZi fa TZi' for i G i#) implies 
c^ij = oti'j for all j G /, although this does not hold for general systems.) 

Since the ordering relationship is transitive and can be extended through- 
out each linkage class by weak reversibility, it follows that the set of reactions 

corresponding to each linkage class Ck either: (1) contributes no stoichio- 

(k) 

metric change to the system (i.e. % G R for all reactions TZi corresponding 
to reactions in C^), or (2) contains a reaction TZi k , G R±, such that 

(k) 

TZj fa TZ ik for all j G B>_ J . We will ignore linkage classes included in case 
(1) since they do not affect (fT2|) . 
We define the set 



k=l 

(fe) 



Now assume that Tjv < 0. This implies that we have 



r 7 v 



> z 



E E w + 

1 ( 

E E + 



it"** 



E 



j&R 



(k) 



Hk - e 



E 



Vi, 



(k) 



(k) 

Regardless of the values of jj, j G R_ , and 7$ fc > for k = 1, . . . 



we 

can pick an e > sufficient small so that Tjv > for every v G K> - In- 
order to satisfy Tjv < for v G M> , therefore, we require v% = for all 

k = !,...,£. We notice that Vi > is permitted 



i G R^ and t G R {k \ 



for i G .Rg^; however, neither Fj nor Tj contain nonzero entries in their 

(k) 

columns corresponding to the elements in i G Rq . It follows that Tv < 
for v G M>o entails v G ker(r). Since this is a sufficient condition for the 
weak dynamical non-emptiability of / by Lemma 13.21 we are done. □ 



13 



In [3], the authors divide semi-locking sets according to whether they 
are critical or not. They first handle the case of non-critical semi- locking 
sets which culminates in Theorem 12.11 It is only in the discussion of critical 
semi-locking sets that they introduce the further condition of dynamical 
nonemptiability. Here we simplify this discussion by showing that every non- 
critical semi-locking set is weakly dynamically non-emptiable and therefore 
falls within the scope of the discussion in Section 13.31 

Theorem 3.4. Consider a semi-locking set I C {l,...,m}. If I is non- 
critical then it is weakly dynamically non-emptiable. 

Proof. Suppose a semi- locking set / C {1, . . . , m} is non-critical. This means 
that / corresponds to the support of a semi-conservation vector so that there 
exists a c € M™ satisfying 

_ ( ci > 0, i e I 

C ~ \ Ci = 0, i<£ I 

such that 

c T r = T . (13) 
It follows from (fl~3|) that there exists a y £ M"q such that 

y T r, = o T . 

Since this implies condition 2. of Farkas' Lemma is satisfied, it follows 
that Condition 1. must necessarily be violated. It follows that any v € M> 
satisfying T/v < must be such that v G ker(r/). By Lemma 13.21 however, 
this is the condition for weak dynamical nonemptiability taking J = which 
is sufficient to prove the result. □ 

3.3 Main Persistence Result 

We are now prepared to present the main result of this paper. The following 
result is a generalization of Theorem 1 2 . 2 1 and includes Theorem l2.1l for mass- 
action kinetics by Theorem 13.41 

Theorem 3.5. Consider a mass-action system with bounded solutions. Sup- 
pose that every semi-locking set is weakly dynamically non-emptiable. Then 
the system is persistent. 
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Proof. We know by Theorem 13.11 that a system with bounded solutions is 
persistent if, for every semi-locking set /, there is an a satisfying 

_ J ai < 0, for i & I 
\ a.i = 0, for i £ I 

so that, for every compact subset K of Lj, there exists a neighbourhood U 
of K in R™ such that (a, f (x)) < for all x G U. 

By assumption, every critical semi-locking set is weakly dynamically non- 
emptiable, which means that there exists an e > and a J C TZj such that 
C(I) n J" e (J) C ker(J, J, e). By LemmaEjS this implies that 

fv < 0, for v G M> v G ker(f). 

It follows that condition 1. of Lemma 13. II is not satisfied. Consequently, in 
order to satisfy 2., there must exist a c G M™Q +nj such that c T f > T . 
We partition c G M> 7 +nj so that 




where c/ G M" 7 and cj G R"^. From this it follows that 

c T f = cjrj + cjrj > T . (14) 

Multiplying through the right-hand side of fjl4[) by i?(x), we have 

cjT/i2(x) + cJr Jj R(x) = -(a, f (x)> + cjrji?(x) > (15) 

where a G R< is the vector — cj extended over the support / and has zeroes 
elsewhere. Clearly a satisfies 

_ J ai < for i G / 
y m = for i g" /. 

By Lemma 12.21 for every compact subset K of Lj and every e > 0, there 
exists a neighbourhood U of in R> such that TjR(x) < 0. It follows 
from (|15p that 

(a,f(x)} <cjrj^(x) <0 

for all x G U. 

Since this holds for every semi-locking set / by assumption, it follows by 
Theorem 13. II that the system is persistent. □ 
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There are several points worth emphasizing about Theorem 13.51 as it 
contrasts with Theorem 12.21 In our result the requirement that the system 
be conservative has been replaced by the more general assumption that 
solutions are bounded, and we do not require the assumption that there 
are no nested critical locking sets. Since a system being conservative implies 
solutions are bounded, the first is not a significant change; however, we have 
opened the result to non-conservative systems for which solutions can be 
bounded by another method, as is the case with complex balanced systems 
(see Section I3.4f) . In Section H] we will see examples where persistence holds 
despite the systems not being conservative. 

We have removed the distinction between critical and non-critical semi- 
locking sets. We do not need to make this distinction since every non-critical 
semi-locking set is weakly dynamically nonemptiable by Theorem 13.41 and 
therefore trivially included in Theorem 13.51 

It is also worth noting that Theorem 13.51 holds even if the onj in ([2]) differ 
from the mass-action exponents in (J3]) provided some regularity conditions 
hold. Specifically, persistence still applies if, rather than (J3j) , we consider 
the reaction vector i?(x) G M> with entries 

m 
3=1 

where Oy > if and only if > and the modified feasibility cone J- e {J) 
is redefined to apply to the partial ordering IZi =3/ TZj with the ct^'s rather 
than the aij. 

3.4 Complex Balanced Systems 

Persistence is of particular interest in the study of the complex balanced 
systems first introduced in 0EJITT] • It is beyond the scope of this paper to 
develop complex balanced systems in any detail except to note that they are 
known to have a unique complex balanced equilibrium concentration within 
each compatibility class and that this concentration is locally asymptotically 
stable relative to the compatibility class (Lemma 4C and Theorem 6A, jllj). 

Despite significant work, however, the following conjecture remains un- 
proven. (This conjecture was first presented in [10]. We state the conjecture 
here as it is stated in [5].) 

Proposition 3.1 (Global Attractor Conjecture). For any complex balanced 
system and any starting point xo € K>o> the associated complex balanced 
equilibrium point x* of C X( , is a global attractor of C xo . 
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Although no general proof is known, many limitations on solutions not 
tending toward x* have been found. Theorem 3.2 of [16] guarantees that 
the w-limit set of a complex balanced system consists either of the unique 
positive equilibrium x* in C Xo or of complex balanced equilibria lying on 
<9R™ ; consequently, persistence of complex balanced systems suffices to af- 
firm Proposition 13.11 Furthermore, since cM> decomposes into the sets 
Lj and w(xo) n Lj ^ for xo G M> implies / is a semi-locking set by 
Lemma [2TTT it follows that we need only prove w(xo) PI Lj = for the sets 
Lj corresponding to semi-locking sets. 

Most recent research on the Global Attractor Conjecture has made use of 
these restrictions ( [TH3"1[5|[T5] ) . For a summary of the major results of this 
research to date, see the discussion preceding Theorem 4.6 of [2]. We will 
append to this result the implications of Theorem 13. 5t however, we begin 
with a Lemma. (In the following we let Fj be associated with the initial 
condition xo G R> an d Fi be associated with xo G R>o-) 

Lemma 3.3. Consider a chemical reaction network. Consider a set I C 
{1, . . . , m} and suppose that Fj ^ for some xo G M>q- Then, for every 
xo G ^>oj cither dim(Fj) = dim(Fj) or Fi = 0. Furthermore, for every 
x G Lj, there exists xq G R> suc ^ that x£F/. 

Proof. Consider the set Lj = {x G M | Xi = if i 6 /}. It is clear that Lj 
is relatively interior to L/, i.e. V x 6 Lj, 3 e > such that -B e (x) n L/ C 
Lj. (B e (x) is the standard Euclidean ball of radius e centered at x.) Now 
consider the affine space (xo + S) n Lj and suppose Fj = (xo + S) n L/ ^ 0. 
Then, VxGF/,3e>0 such that JB £ (x) n [(x + 5) n Lj] = (x + 5) n 
[B e (x) n L/] C (xo + S) (1 Lj = Fj. Consequently, Fj is relatively interior 
to (xo + S) H Lj. Since the dimension of (xo + S) D Lj is the same for all 
xo € M>Q' ft f°ll° ws that dim(F/) is the same for all xo G M> so long as 
Fj 7^ 0. This proves the first claim. 

Since Fj ^ by assumption, we can consider an arbitrary x G i 7 /. 
By definition, we have that (xo — x)j > for i E I and xo — x G 5. 
Now choose an arbitrary x G Lj. It follows from the definition of Lj that 
xo = x + e(xo — x) G KI^q for e > sufficiently small. Since x G Fj for 
xo G M™o an d xGL; was chosen arbitrarily, the second claim follows. □ 

This result guarantees that if Fj is a facet (or vertex) for some xo G 
]R^ j then Fj is a facet (or vertex) for any xo G M> so long as Fi ^ 0. 
Furthermore, it guarantees that L/ can be completely partitioned into sets 
Fj corresponding to facets (or vertices). 
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We are now prepared to prove the following application of Theorem l3.5l to 
complex balanced systems. It should be noted that, while facets are weakly 
dynamically non-emptiable by Theorem 13.31 no comparable result holds for 
vertices (consider the origin in Example 2 of Section U|) . Consequently, the 
following result cannot be attained as a simple application of Theorem 13.51 

Corollary 3.1. Consider a complex balanced mass-action system. Suppose 
that every set Fj corresponding to a semi-locking set I is either a facet, a 
vertex, or empty, or that I is weakly dynamically non-emptiable. Then the 
Global Attractor Conjecture holds for this system. 

Proof. For complex balanced systems, w(xq) D Fj = for every Fj corre- 
sponding to a vertex (Proposition 20 of 0) or the empty set (trivially). Also, 
from Corollary 3.3 of [2], we have that w(xo)nF/ ^ implies u)(xo)ndFj ^ 
for all Fi corresponding to facets; however, since dFi corresponds to some Fj 
not corresponding to a facet, this is a contradiction. It follows from Lemma 
I3.3l that w(xo) PiLj = for any semi-locking set I such that Fj is a facet or a 
vertex for some xo G K>o- It remains to show that w(xo) (~l Lj = for every 
Li corresponding to a weakly dynamically non-emptiable semi-locking set 
I. 

By Theorem [33] we know that for each semi-locking set I which is weakly 
dynamically non-emptiable, there is an a satisfying 



so that, for every compact subset K of Lj, there exists a neighbourhood U 
of K in M™q such that (a,f(x)) < for all x 6 U. Since complex balanced 
systems are bounded, we are justified in using the inductive hypothesis of 
Theorem 3.3 of [15] from \I\ = m to \I\ = 1 to conclude that, for all xo € 
R™ , w(xo) n <9M™ = 0. It follows by Theorem 3.2 of [16] that the Global 
Attractor Conjecture holds for trajectories of such a system, and we are 
done. □ 

4 Examples 

In this section, we present three chemical reaction systems which illustrate 
how Theorem 13.31 Theorem 13.51 and Corollary 13.11 work and what their 
limitations are. 

The first two are examples of non-conservative systems where we exploit 
the results of Section 3. We show that the first example contains a semi- 
locking set which is weakly dynamically non-emptiable but not dynamically 




cti < 0, for i € I 
on = 0, for i I 
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non-emptiable in the sense introduced in [3]- The third is an example of a 
system which does not fall within the bounds of the results discussed in this 
paper or any other papers of which the authors are aware. 

All of the systems considered in this section are complex balanced for all 
sets of rate constants and consequently fall within the scope of the discus- 
sion in Section [3.4i Since the semi- locking set I = {1, . . . , m} corresponds 
to Fj = {0}, which is always a vertex of any compatibility class it is in, we 
may exclude it when considering complex balanced systems, since no tra- 
jectory may approach it. We will call any semi-locking set / = {1, . . . , m} 
trivial. (For further discussion of sufficient conditions to determine complex 
balancing, see [7] and [9].) 



Example 1: Consider the mass-action system 

Ax A 2Ax +A 2 
Ai + A 2 - 

For this system, we have 



r 



o 
-l 



and -R(x) 



k\X\ 
k 2 x\x 2 



and the system is governed by x = TR(x). 

We notice first of all that the system is not conservative and therefore 
does not fall within the scope of the systems considered in [3]. We might 
still be tempted to ask whether the system has semi-locking sets which 
are dynamically non-emptiable, so we consider the semi- locking set I = {1}. 
Relative to this set, we have 1Z 2 1Z\ and 1Z 2 IZ3 so that jF e (I)C\C(I) = 
{0} corresponds to finding a v E K> such that v 1 — v 2 < 0, v 2 < ev\, and 

v 2 < £^3. This can clearly be satisfied for any v = [0 U3] where ^3 > 0. 
Since J- 6 (I) C\C(I) 7^ {0}, it follows that the system contains a critical semi- 
locking set which is not dynamically non-emptiable. 

We notice, however, that Fj is a facet of C Xo = M?,q since s = 2 and 
dim(Fj) = 1. It follows from Theorem 13.31 that I is weakly dynamically 
non-emptiable. Since the system is complex balanced for all sets of rate 
constants and I is the only non-trivial semi-locking set, the Global Attractor 
Conjecture holds for this system by Corollary 13. 11 (This result could also be 
attained by application of Theorem 4.6 of [2], although it should be pointed 
out that Fj is an example of a facet which is not dynamically non-emptiable 
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in the traditional sense so that Corollary 3.5 of the same paper cannot be 
applied.) 



Example 2: Consider the system 



&5 

Ai ^ 2A 2 

fei 

^2 + ^3 & Ai + A 2 . 
The system is governed by the dynamics x = Ti?(x) where 



-1 1-1 1 1 
2-1 0-1 -2 
1-10 



and -R(x) 



k 3 xix 2 

k4X 2 X 3 



This example was first considered in [15], where the authors showed 
that the system is non-conservative, complex balanced for all sets of rate 
constants, and has only the non-trivial semi-locking set / = {1,2}. By the 
methodology presented in that paper, however, they could not find an a 
corresponding to / satisfying ([8]) and ([9]). Since the system is not conser- 
vative, the results of [3] cannot be applied, and since / is not a facet, the 
results of |2] cannot be applied. Here we will show that such an a does in 
fact exist by showing that / is weakly dynamically non-emptiable. 

We have that 



-1 




We have 1Z 2 <i K4, K 3 <i TLi, K 3 <i V.4, and K 5 <i K4 so that 

K T = {(2, 4), (3,1), (3, 4), (5, 4)}. (16) 

We pick the subset J = {(3,4)} so that 

Tj = [ 1 -e ] . 

The condition Tv < for v £ M> is equivalent to the system — v\ + 
v 2 — V3 + V4 + f 5 < 0, 2v\ — v 2 — V4 — 2v 5 < 0, and ^3 — ev^ < for Vi > 0, 
i = 1, . . . , 5. Taking a positive linear combination of these conditions yields 
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v 2 + (1 — 2e)v4 < 0. For < e < 1/2, this can be satisfied for v-i > and 
u 4 > if and only if = 1)4 = 0. It then follows from the third condition 
that V3 = 0. The remaining conditions can be satisfied so long as v\ = v 5 > 
so that 

v€spanj[l 1 ] T | C ker(f ). 

By Lemma [3.21 the semi- locking set / is weakly dynamically non-emptiable. 
Since trajectories are bounded by virtue of the system being complex bal- 
anced, it follows from Theorem 13.51 that the system is persistent and from 
Corollary 13. II that it satisfies Proposition 13.11 

In order to illustrate how the machinery of this result really works, we 
will complete the analysis for / up to the point of applying Theorem 13.11 
From Lemma 13. II we have that there exists ac£ R> such that c T r > T ; 
in fact, we can find it explicitly. This is satisfied if we choose c± = 2, c<i = 1, 
C3 = 2, and < e < 1, for which values we have 



c T f J R(x) =[212 



-1 1 -1 
2-1 
1 



1 1 

-1 -2 
-e 



k\X\ 



-a T x + 2{k^x\X2 — €£14X2X3) > 



where a = [—2 — 1 0] T . It follows that a T ± < 2{k%x\X2 — ££4X2X3) < in 
a neighbourhood of any compact subset of Fj since £3X1X2 < €£4X2X3 under 
the same conditions by Lemma [2.21 This is exactly the condition which was 
expected for application of Theorem 13. 1\ which completes the connection 
with Theorem 13.51 

It is worth reemphasizing that not all sets J satisfying (]16p are sufficient 
to show that / is weakly dynamically non-emptiable. For instance, if we had 
selected J = {(2,4), (3, 1), (3,4)}, we would have had 

1 -e 
-e 1 
1 -e 

In this case, we can satisfy Tv < by choosing 



v G span I [ 1 1 



|T1 



but ker(r) = {0}. Consequently, J is insufficient to show that / is weakly 
dynamically non-emptiable. 
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It is also worth noting that J is not the only choice sufficient for showing 
/ is weakly dynamically non-emptiable. In fact, the maximal set J = TZj 
works with ker(T) = {0}. (In other words, I is dynamically non-emptiable 
in the sense introduced in [3]! We remain unable to use Theorem 4 of [3], 
however, because this system is not conservative.) We can see also that it is 
easier to demonstrate weak dynamical non-emptiability with some choices 
of J than with others, an advantage which would become even more pro- 
nounced for larger systems. 



Example 3: Now consider the system 

Ai + A 2 ^ 3Ai 

ki t ika 

2A 2 <- 2A1 + A3. 
The system is governed by the dynamics x = VR(x) where 



2 -1 
-1 
1 



-2 1 
2 -1 
-1 



and -R(x) 



k\X\X2 

k 2 x\ 
k 3 xfx 3 
k^Xo 



The system is non-conservative, complex balanced for all sets of rate con- 
stants, and has only the non-trivial semi-locking set I = {1, 2}. The system 
is not conservative, so the results of [3] cannot be applied, and / is not 
a facet, so the results of [2] cannot be applied. We consider whether I is 
weakly dynamically non-emptiable. 

We have only the condition 7Z 2 =3/ Ti-z so that J C {(2,3)}. Choosing 
the maximal such set we have 



2 -1 
-1 
1 



-2 1 
2 -1 

-e 



It is clear that v = [0 1 2] satisfies Tv < but 

v ker(f) = span I [-e - e - 1 - 2 + e] T } 

for any e > 0. Since the condition Tv < for v E R> does not imply 
v G ker(r) for the trivial set J = either, it follows that I is not weakly 
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dynamically non-emptiable and thus the results of this paper cannot be 
applied. 

The only other approach that we know of to handle such a situation 
is Corollary 3.15 of [15]. It can be checked, however, that there are eight 
strata which intersect Lj and that the corresponding vectors Ylj=i s /^(j)' 
i = 1, . . . ,5, k = 1, . . . , li — 1, do not have a common a satisfying ([8]) and 
([9]) such that either Condition 1 or Condition 2 is satisfied. We submit, 
therefore, that this is an example of a system whose persistence lies beyond 
the scope of known theory. 
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